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; ABSTRACT 

y—i ■ 

O ' We investigate the cosmological test recently proposed by B. Fort, Y. 

Mellier and M. Dantel-Fort (FMD), where the observed location of the critical 

t— i ■ line in gravitational lensing is used to determine the cosmological parameters, 

^ 1 f2 and A. Applying this method to the cluster of galaxies C10024+1654, FMD 

obtained a constraint on the cosmological constant, A > 0.6, assuming the 

i-C ■ spatially flat universe. It plays a crucial role in this method that the angular 

diameter distance-redshift relation depends on the cosmological models 
through the cosmological parameters. First, using the angular diameter 
distance in the Friedmann-Lemaitre-Robertson- Walker universe, we show 
that one can hardly determine fl by this method without the assumption of 

J^j ■ the spatially flat universe. We also investigate the effect of inhomogeneities 

of the universe by using the Dyer-Roeder angular diameter distance. It is 
shown that the effect of inhomogeneities can become too large to be ignored, 
particularly for a high density universe. As a result, this method cannot 
be taken as a clear cosmological test to determine Q and A, though it may 
provide a bound on Q and A. Moreover, we mention the uncertainty of the 
determination of the velocity dispersion, which is regarded as one of the most 
serious problems in this test. 



Subject headings: cosmology:theory - gravitational lensing, 
galaxies:cluster(C10024+1654) 
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1. Introduction 

In understanding the Universe, it has been one of central issues to determine the 
cosmological parameters; the density parameter Q and the cosmological constant A 
(Weinberg 1972; Peebles 1993). Although some cosmological tests have been proposed 
and applied, the conclusive result has not been obtained yet. The difficulties in 
determining the cosmological parameters consist in the following fact: All cosmological 
tests are based on the fact that relations among observables depend on the cosmological 
parameters. However, the relations also rely on astrophysical models. For instance, the 
galaxy-number count is significantly affected by the formation and evolution of galaxies 
(Yoshii and Takahara 1988; Fukugita, Takahara, Yamashita and Yoshii 1990), and the 
statistics of multiple images of quasars due to the gravitational lensing depends greatly 
on the lens model of galaxies and cluster of galaxies (Flores and Primack 1996; Kochanek 
1996). 

The method to determine the deceleration parameter go, taking type la supernovae 
(SNIa) as the standard candle, has been proposed. Recently, it became possible to apply 
the method to observational data, because of a progress of the search for high redshift 
SNIa (Perlmutter et.al. 1995, 1996; references therein). This method makes use of the 
dependence of the luminosity distance-redshift relation on cosmological parameters. 
It is assumed that the mean absolute magnitude of SNIa is constant over various 
redshifts. The intrinsic dispersion of magnitude of SNIa can be taken into account in 
the statistical sense. Nevertheless, it is likely that the gravitational lensing effect due to 
inhomogeneities of the universe disturbs this test, which has recently been investigated 
by several authors using various approaches (Kantowski, Voughan and Branch 1995; 
Wambsganss et.al. 1996; Frieman 1996). 

On the other hand, the angular diameter distance-redshift relation also depends 
on cosmological parameters. Thus, it is possible to determine cosmological parameters 
through the dependence, if one can measure the redshift and the apparent size of an 
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object whose physical size is known by another independent method. B. Fort, Y. 
Mellier and M. Dantel-Fort have taken notice of the relation between the mass of 
the gravitational lens (cluster of galaxies) and the location of the critical line (Fort, 
Mellier and Dantel-Fort 1996; hereafter FMD). Using the data of cluster of galaxies 
C10024+1654, they obtained a constraint A > 0.6, as a preliminary result. Their 
procedure is as follows (See section 2 for more details): From the observation of a lens 
object (which is nearly axisymmetric along a line of sight) and background galaxies, one 
can determine the velocity dispersion, the angular core radius and the maximum angular 
radius of critical lines. In particular, the critical line, defined as the location where large 
amplification of images occurs (Schneider, Ehlers and Falco 1992, hereafter SEF), can 
be determined observationally as the location where the number density of background 
galaxies is much larger than that in the field (FMD). The location of the critical line 
depends not only on the mass distribution of the lens but also on the ratio between two 
distances; the angular diameter distance from the lens to the source and that from the 
observer to the source. The mass distribution of the lens is usually estimated from the 
velocity dispersion and the core radius. When the redshift of the source is between 3 and 
5, the distance ratio depends strongly on the cosmological parameters but is insensitive 
to the redshift of the source. Thus the cosmological parameters can be determined by 
this dependence. FMD have studied the case of the spatially flat universe. However, it is 
still important to investigate, without assuming the spatially flat universe, how Q and A 
can be really constrained by this method. Moreover, the inhomogeneities of the universe 
may disturb this test by bringing uncertainties to the angular diameter distance-redshift 
relation. We also investigate this effect by using the Dyer-Roeder (DR) angular diameter 
distance (Dyer and Roeder 1972, 1974;hereafter DR). 

This paper is organized as follows. In section 2, the method by FMD is explained. In 
section 3, we consider cases of general Q and A in order to investigate how accurately Q 
and A can be determined by this method. In section 4, the DR angular diameter distance 
is used in order to take into consideration the effect of dumpiness of the universe. In 
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section 5, we discuss how Q and A can be constrained in the clumpy universe. We also 
mention the uncertainty of the determination of the velocity dispersion. Section 6 is 
devoted to summary. 



2. The location of the critical lines and cosmological parameters 

In this section, we would like to summarize the method to determine the cosmological 
parameters by the use of the maximum radius of the critical lines in the gravitational 
lens. It is assumed that the mass profile of the cluster of galaxies is that of the isothermal 
sphere with core (FMD). Then the mass profile of the cluster is written as 



,2 

r 

r 2 /<-. rv2 



2 

p(r) = Po(i + ^)~ , (i) 



where r c is a core radius, po = c /2nGr 2 and a is a line-of-sight velocity dispersion. The 
surface mass density projected to the lens plane is useful in investigating gravitational 
lensing and obtained from Eq.(|l|) as 



a 2 1 



where R is a radius from the center of the cluster on the lens plane. The mass within the 
radius R on the lens plane is expressed as 

M{R) = ^^R* + rl-r c ). (3) 

In the spherically symmetric case such as the isothermal sphere, the critical line appears 
on the Einstein ring (SEF); 



= J™L°%-, (4) 
V c 2 D 0L D s 

where 9 cr and 9e are the angular radius of the critical line and that of the Einstein ring, 
and Dol, Dos and Dls are the angular diameter distances between the observer and 
the lens, between the observer and the source, and between the lens and the source, 
respectively. In the case of isothermal sphere, we may use Eq.@ to obtain 



V c 2 v ' DqlDos 
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where R cr is the radius of the critical line. It should be noted that the observable is not 
the radius of the critical line but its angular radius, 9 cr = R cr jD OL . From Eq.(^[), we 
obtain the important ratio (FMD), which depends only on the observables, 

D LS C 2 



((e+a i/2 +^), (6) 



D os 47ra 2 

where we defined 8 C = t c /Dql- Usually, the radius of the critical line (several tens arcsec, 
several Mpc) is much larger than the core radius (less than Mpc). Then we can use the 
approximate relation 

° LS - * i„. (7) 



Dos 47TCT 2 

It is noteworthy that the angular diameter distance has not been specified. Since the 
angular radius of the critical line 8 cr , the line-of-sight velocity dispersion in the cluster a 
and the angular core radius 9 C can be obtained observationally, the ratio Dls/Dos can 
be determined from these data. On the other hand, the distance ratio Dls/Dos as we U 
as the angular diameter distances, D LS and D Q s, depends on the cosmological models. 
For instance, in the FLRW universe, the relation between the redshift and the distance 
ratio depends only on the cosmological parameters, Q and A (Weinberg 1972; Peebles 
1993). Figure 1 shows how the relation between the distance ratio Dls/Dos and the 
redshift of the source depends on the cosmological parameters, where we assume the 
angular diameter distance in the FLRW universe. It also shows that Dls/Dos increases 
with the redshift of the source. It should be noted in Fig.l that D LS /D s is sensitive 
to Q and A, but insensitive to the redshift of the source from 3 to 5. Therefore, if one 
obtains observationally the maximum of D is /Dos, the cosmological parameters can be 
constrained without the precise measurement of redshift of background faint galaxies 
(FMD). In order to constrain the cosmological parameters, it is necessary to obtain 
Dls/Dos within about ten percent accuracy, for C10024+1654 at zl = 0.39, as shown 
by Fig.l. 

From Eq.(£|), it is clear that the cosmological test discussed here relies on (1) which 
angular diameter distance to be used and (2) the mass of the cluster as the gravitational 
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lens. We mainly investigate the effect of the first aspect on the cosmological test. In 
particular, we pay much attention to the distance ratio D LS /D s- It is worthwhile 
to mention that the cosmological test discussed here does not depend on the Hubble 
constant because of our using only the distance ratio Dls/Dos- 

3. Cosmological test in the FLRW universe 

In the previous section, we have given Fig.l, which shows a relation among zs, 
D LS /D os and cosmological parameters, by assuming Q + A = 1. In the remaining of this 
paper, we fix the redshift of the source at z s = 4 in order to make our discussion clear. 
This treatment is plausible in practice, since D LS /D s is insensitive to the redshift of the 
source for zs > 3, as shown by Fig.l. Figure 2 shows relations between Q and Dls/Dos 
in two cases of Q + A = 1 and A = 0. It is found that a bound on A can be obtained in 
the case of f2 + A = 1, particularly for large D LS /D s- Remarkably, D LS /D s depends 
very weakly on Q when A = 0. However, it is not clear whether Dls/Dos is generally 
insensitive to Q. Thus we investigate the dependence of Dls/Dos on Q and A. In order 
to investigate a relation among Q, A and D L s/D s, we plot contours of D L s/D s on 
fi-A plane assuming zl = 0.5, as shown by Fig. 3. Figure 3 shows clearly that D L s/D s 
is very insensitive to Q regardless of A. Therefore, it is difficult to determine Q by this 
method. However, when we assume the spatially flat universe, one can put a bound on 
the cosmological constant, as FMD discussed. We use Fig. 3, in order to illustrate how 
we can constrain the cosmological constant for a given Dls/Dos assuming Q + A = 1. 
For instance, Dls/Dos ~ 0.80 would suggest a universe dominated by the cosmological 
constant (A ~ 0.9). On the other hand, if Dls/Dos is less than 0.67, the vanishing 
cosmological constant is preferred. 

Figs.4(a)-(c) show contours of Dls/Dos on the Q-\ plane in the cases of zl = 0.1, 0.3 
and 1.0 respectively . These suggest that, regardless of the redshift of the lens, D L s/D s 
depends very weakly on Q, but is more sensitive to A in a lower density universe. Thus, 
one may constrain the cosmological constant by using the lens at various redshifts 
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as well. However, as shown by Fig.4(a), the cosmological constant does not change 
D L s/Dos so much for a low redshift lens, that one must use a lens at a moderately 
high redshift (zl > 0.3) for this cosmological test. It is important to note the difference 
between the new method and the so-called qo test; the new method is more sensitive to 
A than to Q, while the q test aims at the determination of the combination VL/2 — A in 
the matter-dominated universe (Weinberg 1972; Peebles 1993). 



On a sufficiently large scale, our universe can be described by the FLRW universe. 
Nevertheless, the light ray propagates from the source to the observer according to the 
geodesic equation which is determined by the local metric along the light ray (Weinberg 
1972; Peebles 1993). Therefore, it is likely that the light propagation is affected by 
the inhomogeneities of the universe. Much attention has focused on the effect of the 
inhomogeneities on the luminosity or angular diameter distance between the observer 
and the source, Dos (DR; Futamase and Sasaki 1989; Kasai, Futamase and Takahara 
1990; SEF; Sasaki 1993). Thus, here we must consider how much the inhomogeneities 
affect the ratio of angular diameter distances D LS /D s- 

We consider the DR angular diameter distance, for its simplicity, in order to take 
into account the inhomogeneities (DR; SEF). The DR angular diameter distance is 
determined by 



where the parameter a represents the dumpiness along the light ray. In the FLRW 
universe, a is unity, while a vanishes when the light ray propagates through the empty 
space. In this paper, we consider < a < 1, since we have little knowledge of a. Here 
w is an affine parameter, which is assumed (SEF) to be given by that in the FLRW 
universe 



4. Cosmological test in the clumpy universe 




dw 2 



(8) 
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It is noteworthy that the coefficient of the last term of Eq.(P), 3afi/2, comes from the 
Ricci focusing by the clumpy matter and, therefore, the DR angular diameter distance is 
a decreasing function of a for a fixed redshift (DR; SEF). 

Since we have considered the FLRW universe (a = 1) in the previous section, here 
we consider the other limiting case a = in order to investigate the effect of dumpiness 
of the universe. 

4.1. The dependence on and A 

By solving Eq.([8]) numerically, we can obtain the angular diameter distances Dls 
and Dos- As a result, we draw Fig.5 which shows relations between Q and Dls/Dqs in 
two cases of Q + A = 1 and A = 0. It is found that, in the clumpy universe, the distance 
ratio Dls/Dos is smaller than that in the case of a = 1. In particular for a higher 
density universe, the dependence of Dls/Dos on a is larger in both cases; Q + A = 1 
and A = 0. Next, we draw contours of Dls/Dqs on Q-X plane in Fig.6, in order to 
clarify the dependence of Dls/Dos on Q and A in the clumpy universe. By comparing 
Fig.6 with Figs. 3 and 4, we find that, in the universe with any Q and A, Dls/Dqs for 
a = is smaller than that for a = 1. For any A, we also find that the discrepancy of 
Dls/Dqs between a = and 1 increases with Q. This may be attributed to the fact 
that the focusing effect due to the dumpiness in Eq.(||) increases with Q. Furthermore, 
Fig.6 implies that, for a = 0, Dls/Dqs is very sensitive to Q rather than to A. This 
is in contrast to the case of a = 1. We will show an important property of Dls/Dos 
in the next subsection, before we will discuss in section 5 how one can constrain the 
cosmological parameters by taking into account the effect due to the clumpy universe. 
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4.2. The general case 

In the above, we have only considered two limiting cases of a = and 1. In Figs. 2 
and 5, the following interesting property of the distance ratio appears 

^-(a = 0; Empty) < ^-(a = 1; FLRW). (10) 
Dos D s 

By noticing that a enhances monotonically the Ricci focusing in Eq.@, one expects 
that Eq. (|T0|) is generalized to a relation for a\ < a<i as 

Dls , x ^ D LS . . , , 

TJ— ("i) < 7T-(«2). (11) 

u OS u OS 

This is proved as follows: Let us fix the redshift of the source and the cosmological 
parameters, Q and A. Then we define a function of zl parameterized by a as 

X a (z L ) = *£*(a). (12) 
Uos 

From Eq.(§), we obtain the equation for X a (zi) 

-£rX a (z L ) + |(i + ^) 5 «fix Q (z L ) = 0, (13) 

where u>£ is an affine parameter at the lens. We define the Wronskian as 

W(X„,X m ) = {x„^-xJ-^). (14) 
Using Eq. (fD|) , we obtain 

J-W(X ai ,X a2 ) <0, (15) 
where we used ai < a 2 . Since both X ai and X„ 2 vanish at zl = z s , we obtain 

W(X ai (z s ),X a2 (z s )) = 0. (16) 
From Eqs.(|l5]) and (fEBI), we find 

W(X ai ,X aa )>0, (17) 
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where we used the fact that the affine parameter w defined by Eq.([J) is an increasing 
function of z. Eq.(|I7|) is rewritten as 

-^ln^>0. (18) 
dw L X ai 

Since X a always becomes unity at the observer, we find 

fa *-f t = °> =0. (19) 
X ai (z L = 0) 



ln^>0. (20) 



From Eqs.(|T8D and (H), we find 



This gives us what we wish to prove. It should be noted that Eq.(]20|) holds even if one 
uses the opposite sign in the definition of the affine parameter in Eq.(^). Eq.(|ll|) means 
that the distance ratio D LS /D s for < a < 1 takes a value between those for a = 
and 1, namely for < a < 1, 

^(o = 0; Empty) < ^£(«) < ^(a = 1; FLRW) . (21) 
J^os u OS ^os 

Therefore, all we have to do is to investigate two limiting cases a = and 1, in order 
to know a region where Dls/ Dos for < a < 1 takes its value for given zl, Q and A. 
Eq.(PTD shows that, for the light propagation through a lower density region such as a 
void, the distance ratio D LS /D s is always less than that in the FLRW universe. 



5. Discussion 



We have used the DR angular diameter distance, in order to investigate how much 
the clumpy universe affects the determination of the cosmo logical parameters by the 
method using the location of the critical lines in the gravitational lens. It is shown in 
section 4 that the distance ratio Dls/ Dos m a clumpy universe is necessarily less than 
that in the FLRW universe. It is also shown that the distance ratio Dls I Dos m a higher 
density universe is more affected by the dumpiness. As a result, it seems impossible 
that the cosmological parameters can be clearly determined by this method. This 
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situation may remind us of the recent discussion about the q test using the luminosity 
distance-redshift relation, where about fifty-percent deviation can happen in the Swiss 
cheese model of the universe (Kantowski, Vaughan and Branch 1995). 

Nevertheless, there are some cases where the cosmological parameters can be 
constrained by noticing that the dumpiness always decreases Dls/Dos- In order to 
illustrate these cases, we use Figs. 3 and 6(c) [z L = 0.5), assuming two cosmological 
models with (1) Q + A = 1 and (2) A = 0. First, we consider the case of (1) Q + A = 1: 
In the case of Dls/Dos between 0.67 and 0.85, one can put a bound on A, though the 
bound is sensitive to a. The allowed region of A becomes wide as Dls/Dos decreases. 
If Dls/Dos is between 0.53 and 0.67, an upper bound on A can be obtained, which is 
given by a = 0. Next, we use the case of (2) A = 0: In the case of Dls/Dos > 0.63, 
one can put a lower bound on Q, which is given by a = 1. It is noteworthy that Q must 
be greater than unity if Dls/Dos exceeds 0.67. If Dls/Dos is less than 0.63, one can 
also put a lower bound on Q, which is given by a = 0. In the above, we assumed either 
fi + A = lorA = 0a priori. However, even without such a constraint, the allowed region 
of the parameters Q and A is still bounded by two limiting cases; a = and 1 (For 
instance, see two curves for Dls/Dos — 0.69 in Fig.6(c)). It is important to mention 
that a lower bound on A, given by using the angular diameter distance in the FLRW 
universe (a = 1), is not affected by the dumpiness. Thus, A > 0.6 given by FMD cannot 
be relaxed by taking into account the dumpiness of the universe. Here, it should be 
noted that the cosmological constant may be well constrained by this method, if the 
density parameter can be determined by another method. 

In the above, we have considered only the dependence of Dls/Dos on the 
cosmological parameters. Here, we briefly mention the dependence of this test on 
astrophysical models of the cluster of galaxies. As shown in Eq.(^|), the precise 
determination of the velocity dispersion is necessary even if one assumes the isothermal 
sphere as the lens. If one specifies the redshifts of the lens and the source, the accuracy 
required for the determination of the velocity dispersion depends on which cosmological 
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models we wish to compare with. If one wishes to distinguish, for instance, the 
cosmological model with (O, A) = (0.2,0) from that with (ft, A) = (0.2,0.8), one must 
distinguish Dls/Dos = 0.64 from 0.76 for zl = 0.5 and zs = 4. This means that one 
must discriminate the difference of Dls/Dos at worst by 16 percent. However, the 
statistical error of the velocity dispersion is at least about ten percent because of the 
limited number(about a hundred) of galaxies inside the radius of the critical line, even 
though the velocity of individual galaxies can be measured precisely (Dressier, Gunn 
and Schneider 1985). This required accuracy is almost unchanged around z L = 0.5. 
Thus, because of the statistical error of the determination of the velocity dispersion, it 
is concluded that the current data of C10024+1654 does not rule out the cosmological 
model without A. However, in the statistical sense, it is possible to put a bound on ft 
and A, if one can apply this method to a lot of clusters of galaxies. For instance, if 
one uses ten clusters of galaxies, the statistical error of the velocity dispersion becomes 
less than four percent, so that one may distinguish some cosmological models, say with 
(ft, A) = (0.2,0) and (0.2,0.8). 

Furthermore, the mass profile used as the lens model must be confirmed in order to 
make this cosmological test reliable. Thus it is of great importance that the mass or the 
mass profile of the lens is determined precisely. It is hoped that the mass reconstruction 
by the weak gravitational lensing (Kaiser and Squires 1993; Kaiser 1995; Seitz and 
Schneider 1995a, 1995b) as well as the analysis of the velocity field in the cluster of 
galaxies (Yee, Ellingson and Carlberg 1996) clarifies the precise mass profile of the cluster 
as the lens. Then the cosmological test discussed here will work well. 

6. Summary 

Without assuming the spatially flat universe, we have reexamined the cosmological 
test done by FMD, by using the DR angular diameter distance in order to take into 
account the dumpiness of the universe. If the angular diameter distance is that in 
the FLRW universe, namely ignoring the effect of the dumpiness, it is found that ft 
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cannot be determined directly by this method. It was also shown that the effect of 
inhomogeneities of the universe can become too large to be ignored, particularly for a 
high density universe. Therefore, this method cannot be taken as a clear cosmological 
test to determine Q and A. However, it is possible to provide a bound on Q and A in 
the clumpy universe. For more quantitative discussion, it is necessary to consider the 
gravitational lensing in a more realistic inhomogeneous universe. 

In section 5, we mentioned that the statistical error of the velocity dispersion can be 
the most serious problem in the cosmological test discussed here. Therefore, in order to 
constrain Q and A by this method, it is necessary to observe many clusters of galaxies. 
For instance, Sloan Digital Sky Survey (SDSS) will provide us with a lot of clusters of 
galaxies which are useful for this cosmological test. 

The author would like to thank M. Sasaki and T. Tanaka for careful reading of the 
manuscript and invaluable comments on it. He is also very grateful to B. F. Roukema 
for pointing out the recent work by B. Fort, Y. Mellier and M. Dantel-Fort. He would 
like to thank S. Ikeuchi and M. Sasaki for useful advice and continuous encouragement. 
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Fig. 1— The distance ratio D LS /D os for the lens at z L = 0.4. (SI, A) = (1,0), 
(f2, A) = (0.4, 0.6) and (Q, A) = (0, 1) are shown by solid, dashed and dotted lines 
respectively. 

Fig. 2. — The distance ratio D LS /D os for the lens at z L = 0.5 and the source at z s = 4. 
We use the angular diameter distance in the FLRW universe. We assume f2 + A = 1 (solid 
line) and A = (dashed line). 

Fig. 3. — Contours of D LS /D os on Q-X plane Each contour from top to bottom is drawn 
where Dls/Dos runs from 0.85 to 0.65 at the interval of 0.02. The dotted line connecting 
(O, A) = (1, 0) and (0, 1) indicates the spatially flat universe. 

Fig. 4. — Contours of Dls/Dos on plane using the angular diameter distance of 
the FLRW universe. The dotted line connecting (Q, A) = (1, 0) and (0, 1) indicates the 
spatially flat universe, (a) We assume zl = 0.1 and zs = 4. Each contour from top to 
bottom is drawn where D LS /D Q s runs from 0.97 to 0.91 at the interval of 0.01. (b)We 
assume zl = 0.3 and zs = 4. Each contour from top to bottom is drawn where Dls/Dos 
runs from 0.90 to 0.76 at the interval of 0.02. (c)We assume zl = 1 and z s = 4. Each 
contour from top to bottom is drawn where D L s/D s runs from 0.64 to 0.44 at the 
interval of 0.02. 

Fig. 5. — The distance ratio D L s/D s for the lens at zl = 0.5 and the source at zs = 4. 
We use DR angular diameter distance with a = 0. We assume f2 + A = 1 (solid line) and 
A = (dashed line). 

Fig. 6. — Contours of Dls/Dos on Q-X plane, using DR angular diameter distance with 
a = 0. The dotted line connecting (f2, A) = (1,0) and (0, 1) indicates the spatially flat 
universe, (a) We assume z L = 0.1 and z s = 4. Each contour from top to bottom is drawn 
where D L s/D s runs from 0.97 to 0.87 at the interval of 0.01. (b)We assume zl = 0.3 
and zs = 4. Each contour from top to bottom is drawn where Dls/Dos runs from 0.90 
to 0.68 at the interval of 0.02. (c)We assume z L = 0.5 and z s = 4. Each contour from top 
to bottom is drawn where Dls/Dqs runs from 0.85 to 0.53 at the interval of 0.02. The 
dashed line is drawn at the level of Dls/Dos — 0.69 for a = 1 for comparison. (d)We 
assume zl — 1 and z$ = 4. Each contour from top to bottom is drawn where Dls/Dqs 
runs from 0.70 to 0.34 at the interval of 0.02. 
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